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NUMERICAL STUDIES OF WEAKLY STOCHASTIC MAGNETIC 

RECONNECTION 

G. Kowal,^'^ A. Lazarian/ E. T. Vishniac^ and K. Otmianowska-Mazur^ 

RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to translate your abstract into 
QQ \ Spanish, the editors will do it for you. We study the effects of turbulence on magnetic reconnection using 

■ three-dimensional numerical simulations. This is the first attempt to test a model of fast magnetic reconnection 
\ proposed by Lazarian & Vishniac (1999), which assumes the presence of weak, small-scale magnetic field 

■ structure near the current sheet. This affects the rate of reconnection by reducing the transverse scale for 
O I reconnection flows and by allowing many independent flux reconnection events to occur simultaneously. We 

■ performed a number of simulations to test the dependencies of the reconnection speed, defined as the ratio 
\ of the inflow velocity to the Alfven speed, on the turbulence power, the injection scale and resistivity. Our 

^-H ■ results show that turbulence significantly affects the topology of magnetic field near the diffusion region and 
\ increases the thickness of the outflow region. We confirm the predictions of the Lazarian & Vishniac model. 

■ In particular, we report the growth of the reconnection speed proportional to ~ F^^, where Vi is the amplitude 
\ of velocity at the injection scale. It depends on the injection scale linj as ~ {linj / L)'^^^ ^ where L is the size of 

Qh' the system, which is somewhat faster but still roughly consistent with the theoretical expectations. We also 
Q \ show that for 3D reconnection the Ohmic resistivity is important in the local reconnection events only, and the 
^ ■ global reconnection rate in the presence of turbulence does not depend on it. 

: ABSTRACT 

We study the effects of turbulence on magnetic reconnection using three-dimensional numerical simulations. 

■ This is the first attempt to test a model of fast magnetic reconnection proposed by Lazarian & Vishniac 
^ ! (1999), which assumes the presence of weak, small-scale magnetic field structure near the current sheet. This 

affects the rate of reconnection by reducing the transverse scale for reconnection flows and by allowing many 
', independent flux reconnection events to occur simultaneously. We performed a number of simulations to test 
(S| ■ the dependencies of the reconnection speed, defined as the ratio of the inflow velocity to the Alfven speed, 
on the turbulence power, the injection scale and resistivity. Our results show that turbulence significantly 
affects the topology of magnetic field near the diffusion region and increases the thickness of the outflow region. 
OO ! We confirm the predictions of the Lazarian & Vishniac model. In particular, we report the growth of the 

■ reconnection speed proportional to ~ F^^, where Vi is the amplitude of velocity at the injection scale. It 
^ \ depends on the injection scale linj as ~ (linj / L^^^ ^ where L is the size of the system, which is somewhat 

■ faster but still roughly consistent with the theoretical expectations. We also show that for 3D reconnection the 
rS \ Ohmic resistivity is important in the local reconnection events only, and the global reconnection rate in the 

^ ■ presence of turbulence does not depend on it. 
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1. INTRODUCTION 

Magnetic fields play a key role in the astrophys- 

ical processes such as star formation, the transport 
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Reconnection is a fundamental process that de- 
scribes how bundles of magnetic field lines can pass 
through each other. Different approaches to the prob- 
lem of recon nection are discussed in th e compan- 
ion paper by iLazarian fc Vishniad ( 2QQ8I ). In what 
follows we concentrate on testing the model of re- 
connection for a weakly sto chastic field proposed by 
Lazarian fc Vishniac (|l999l . LV99 henceforth). They 
argued that reconnection speed is equal to the up- 
per limit imposed by large-scale field line diffusion, 
expressed by 



Vrec = Va min 



1/2 



1/2- 



Yl 

Va 



1/2 



(1) 



where Va is the Alfven speed, L is the size of the 
system, / is the injection scale, and Vi is the velocity 
amplitude at the injection scale. In this relation, the 
reconnection speed is determined by the characteris- 
tics of turbulence, namely, its strength and injection 
scale. Most importantly, there is no explicit depen- 
dence on the Ohmic resistivity. 

The numerical testing of reconnection models is 
far f rom trivial. While mos t of the reconnection work 
(see Priest fc ForbesI |2QQQ[ ) is performed in 2D, the 
LV99 model is intrinsically 3D. 

In ^we describe our numerical model, in ^we 
describe the reconnection rate we use, in §31 and ^Sl 
we describe and discuss the results we obtained, and 
in §6] we summarize the paper drawing the main con- 
clusions. 

2. NUMERICAL MODELING OF LV99 
RECONNECTION 

2.1. Governing Equations 
We use a higher-order shock-capturing 
Godunov-type scheme based on the essentially 
non oscillatory (ENO) spacial reconstruction 
and Run^e-Kutta (RK) time integra tion (see 
Del Zanna. Bucciantini fc Londrillol 2QQ3I , e.g.) to 
solve isothermal non-ideal MHD equations, 

dp 

dt 



T^ + V-M = 0, 



dpv 



pvv + ptI 



BB 

in 



OA 

'dt 



E 



0, 



(2) 
(3) 
(4) 



where p and v are plasma density and velocity, re- 
spectively, A is vector potential, E = —vx B-\-r]j is 
electric field, ^ = V x ^ is magnetic field, j =\/ x B 
is current density, pr = a'^ p-\- B'^ / Sir is the total pres- 
sure, a is the isothermal speed of sound, r] is resis- 
tivity coefficient, and / represents the forcing term. 



We incorporated the field int erpolated c on- 
strained transport (CT) scheme (see Toth 2QQQ[ ) in 
to the integration of the induction equation to main- 
tain the V • -B = constraint numerically. 

Some selected simulations that we perform in- 
clude anomalous resistivity modeled as 

V = Vu + vJ^-i)h(^), (5) 

VJcrit / VJcrit / 

where r]u and r]a describe uniform and anomalous re- 
sistivity coefficients, respectively, jcrit is the critical 
level of the absolute value of current density j above 
which the anomalous effects start to work, and H is 
a step function. For most of our simulations r]a = 0^ 
however. 

2.2. Initial Conditions and Parameters 

Our initial magnetic field is a Harris current sheet 
of the form Bx = Bxq tanh(?//6>) initialized using the 
magnetic vector potential Az = \n\cosh{y/0)\. In 
addition, we use a uniform shear component Bz = 
Bzo = const. The initial setup is completed by set- 
ting the density profile from the condition of the uni- 
form total pressure pt(^ = 0) = const and setting the 
initial velocity to zero everywhere. 

In order to initiate the magnetic reconnection we 
add a small initial perturbation of vector potential 
SAz = Bxocos{27rx) ex.-p[—{y/d)'^] to the initial con- 
figuration of Az{t = 0). The parameter d describes 
the thickness of the perturbed region. 

Numerical model of the LV99 reconnection is 
evolved in a box with open boundary conditions 
which we describe in the next sub-section. The box 
has sizes Lx = Lz = 1 and Ly = 2 with the res- 
olution 256x512x256. It is extended in Y-direction 
in order to move the infiow boundaries far from the 
injection region. This minimizes the infiuence of the 
injected turbulence on the infiow. 

Initially, we set the strength of anti-parallel mag- 
netic field component to 1 and we vary the shear 
component Bz between 0.0 and 1.0. The speed of 
sound is set to 4. In order to study the resistiv- 
ity dependence on the reconnection we vary the re- 
sistivity coefficient r] between values 0.5 • 10~^ and 
2 • 10~^ which are expressed in dimensionless units. 
This means that the velocity is expressed in units 
of Alfven speed and time in units of Alfven time 
tA = L/Va^ where L is the size of the box. 

2.3. Boundary Conditions 
The boundary conditions are set for the ffuid 
quantities and the magnetic vector potential sepa- 
rately. For density and velocity we solve a wave equa- 
tion with the speed of propagation equal to the max- 
imum linear speed which is the fast magnetosonic 
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one and its sign corresponding to the outgoing wave. 
This assumption guarantees that ah waves generated 
in the system are free to leave the box without sig- 
nificant reflections. Moreover, the waves propagate 
through the boundary with the maximum speed, 
reducing the chance for interaction with incoming 
waves. 

For the vector potential, its perpendicular com- 
ponents to the boundary are obtained using the first 
order extrapolation, while the normal component 
has zero normal derivative. This guarantees that the 
normal derivative of the magnetic field components 
is zero, which reduces the influence of the magnetic 
field at the boundary on the plasma flow. 

This type of boundary conditions represents a 
mixed inflow/outflow boundaries, which are adjust- 
ing during the evolution of the system. It means that 
we do not set fixed values of quantities and do not 
drive the flow at the boundaries in order to achieve 
a stationary reconnection. 

2.4. Method of Driving Turbulence 
In our mo del we dr i ve tur bulence using a method 
described by lAlvehusI (|l999h . The forcing is imple- 
mented in spectral space where it is concentrated 
with a Gaussian profile around a wave vector cor- 
responding to the injection scale linj. Since we can 
control the scale of injection, the energy input is in- 
troduced into the flow at arbitrary scale. The ran- 
domness in the time makes the force neutral in the 
sense that it does not directly correlate with any of 
the time scales of the turbulent flow and it also makes 
the power input determined solely by the force-force 
correlation. This means that it is possible to gen- 
erate different states of turbulence, such as axisym- 
metric turbulence, where the degree of anisotropy 
of the forcing can be chosen a priori through the 
forcing parameters. In the present paper we limit 
our studies to the isotropic forcing only. The total 
amount of power input from the forcing can be set 
to balance a desired dissipation at a statistically sta- 
tionary state. In order to get the contribution to 
the input power in the discrete equations from the 
force-force correlation only, the force is determined 
so that the velocity-force correlations vanish for each 
Fourier mode. 

On the right hand side of Equation ([3]), the forc- 
ing is represented by a function / = pa, where p 
is the local density and a is a random acceleration 
calculated using the method described above. 

The driving is completely solenoidal, which 
means that it does not produce density fluctuations. 
Density fluctuations results from the wave interac- 
tion generated during the evolution of the system. 
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Fig. 1. Evolution of total mass M and kinetic and mag- 
netic energies, Ekin and Emag respectively. The kinetic 
energy Ekin has been amplified by a factor of 100 to vi- 
sualize its evolution more clearly. The resistivity in this 
model is ?7 = 10""^ and the shear component of magnetic 
field = 0.1. 

Nevertheless, in our models we set large values of 
the speed of sound approaching nearly incompress- 
ible regime of turbulence. 

3. RECONNECTION RATE MEASURE 

We measure the reconnection rate by averaging 
the inflow velocity Vin divided by the Alfven speed 
Va over the inflow boundaries. In this way our defi- 
nition of the reconnection rate is 



Vrec = {Vin/VA)s 



y — yminiymax 



V 



dS, (6) 



where S defines the XZ planes of the inflow bound- 
aries. 

4. RESULTS 

In this section we describe the results obtained 
from the three dimensional simulations of the mag- 
netic reconnection in the presence of turbulence. 
First, we investigate the Sweet-Parker reconnection, 
a stage before we inject turbulence. A full under- 
standing of this stage is required before we perform 
further analysis of reconnection in the presence of 
turbulence. 

4.1. Sweet-Parker Reconnection 



As we described in §2.2[ Sweet-Parker reconnec- 
tion develops in our models as a result of an initial 
vector potential perturbation. In order to reliably 
study the influence of turbulence on the evolution of 
such systems, we need to reach the stationary Sweet- 
Parker reconnection before we start injecting turbu- 
lence. 

Figure [1] shows the evolution of total mass, ki- 
netic and magnetic energies until the moment at 
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which we start injecting the turbulence (i.e. t = 7). 
All shown quantities, after some initial adaptation, 
reach steady, almost constant values. The near zero 
time derivatives of total mass, kinetic and magnetic 
energies guarantee the stationarity of the system. 
We remind the reader, that the system evolves in the 
presence of open boundary conditions, which could 
violate the conservation of mass and total energy. 
As we see, the conservation of these quantities is 
well satisfied during the Sweet-Parker stage in our 
models. 



0.10 7 
0.08 - 




Time [Ta = L / Vj 

Fig. 2. Evolution of the reconnection rate (Vin/VA) for 
the same model as in Fig.[T] The reconnection rate grows 
initially until it reaches the stationary solution. 

The reconnection rate, shown in Figure [21 also 
confirms that the evolution reaches a stationary 
state. Initially, the reconnection rate {Vin/VA) grows 
until time t ~ 1.8, when it reaches the maximum 
value of ~ 0.06. Later on, it drops a bit approach- 
ing a value of 0.04. During the last period of about 
3 Alfven time units, the change of the reconnection 
rate is very small. We assume that these conditions 
guarantee a nearly steady state evolution of the sys- 
tem, so at this point we are ready to introduce tur- 
bulence. 

In Figure [3] we present the velocity and magnetic 
field configuration of the steady state. In the left 
panel of the figure we show the topology of veloc- 
ity field as textures. The brightness of texture cor- 
responds to the amplitude of velocity. The texture 
itself shows the direction of the field lines. The topol- 
ogy of the velocity field is mainly characterized by 
strong outflow regions along the mid-plane. The out- 
flow is produced by the constant reconnection pro- 
cess at the diffusion region near the center and the 
ejection of the reconnected magnetic flux through 
the left and right boundaries. The system is in a 
steady state when the flux which reconnects is coun- 
terbalanced by the incoming fresh flux. The inflow is 
much slower then the outflow, but still its direction 



can be recognized from the texture shown in the left 
plot of Figure [3l 

The topology of the magnetic field is presented 
in the middle panel of Figure [3l We recognize the 
anti-parallel configuration of the field lines with the 
uniform strength out of the mid-plane. Near the 
mid-plane, the horizontal magnetic lines are recon- 
nected generating the Y-component, which is ejected 
by the strong outflow. In addition, we show the abso- 
lute value of current density in the right panel of Fig- 
ure [3l We see very elongated diffusion region in the 
middle of the box, where reconnection takes place. 
The maximum value of \j\ does not exceed a value 
of 25. 

4.2. Effects of Turbulence 

The essential part of our studies covers the ef- 
fects of turbulence on reconnection. Our goal is 
to achieve a stationary state of Sweet-Parker recon- 
nection, which is described in the previous subsec- 
tion, and then introduce turbulence at a given injec- 
tion scale linj, gradually increasing its strength to 
the desired amplitude corresponding to the turbu- 
lent power Pinj' We inject turbulence in the region 
surrounding the mid plane and extending to the dis- 
tance of around one quarter of the size of the box. 
The transition period during which we increase the 
strength of turbulence, has length of one Alfvenic 
time and starts at t = 7. It means that from t = 8 
we inject turbulence at maximum power Pinj- 

In Figure m we show examples of XY-cuts (upper 
row) and XZ-cuts (lower row) through the box of 
the velocity (left panel) and magnetic field (middle 
panel) topologies with the intensities corresponding 
respectively to the amplitude of perpendicular com- 
ponents of velocity and magnetic field to the normal 
vector defining the plotted plane. 

The first noticeable difference compared to the 
Sweet-Parker configuration is a significant change of 
the velocity and magnetic field topologies. Velocity 
has very complex and mixed structure near the mid 
plane, since we constantly inject turbulence in this 
region (see the left panel in Fig. Hj). Although the 
structure is very complex here, most of the velocity 
fluctuations are pointed in the directions perpendic- 
ular to the mean magnetic field. This comes from 
the fact, that in the nearly incompressible regime 
of turbulence, most of the fluctuations propagate as 
Alfven waves along the mean magnetic field. Slow 
and fast waves, whose strengths are significantly re- 
duced, are allowed to propagate also in directions 
perpendicular to the mean field. As a result most of 
the turbulent kinetic energy leaves the box along the 
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Fig. 3. Topology and strength of velocity field (left panel) and magnetic field (middle panel) during the Sweet-Parker 
reconnection at t = 7. In the right panel we show the absolute value of current density j. The images show the XY-cut 
through the domain at Z = 0. 
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Fig. 4. Topology and strength of velocity field (left panel) and magnetic field (middle panel) in the presence of fully 
developed turbulence at time t = 12. In the right panel we show distribution of the absolute value of current density | J|. 
The images show the XY-cut (upper row) and XZ-cut (lower row) of the domain at the midplane of the computational 
box. Turbulence is injected with the power Pinj = 1 at scale kinj = 8. 
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Fig. 5. Time evolution of the reconnection speed Vrec 
for models with different powers of turbulence Pinj (see 
legend) . In the subplot we show the dependence of Vrec 

on Pinj. 

magnetic lines. The fluctuations, however, efficiently 
bend magnetic lines near the diffusion region. There, 
the strength of the magnetic field is reduced, since 
this is the place where magnetic lines change their di- 
rections (see the middle upper panel in Fig. [4]). The 
interface between positively and negatively directed 
magnetic lines is much more complex then in the 
case of Sweet-Parker reconnection. This complexity 
favors creation of enhanced current density regions, 
where the local reconnection works faster since the 
current density reaches higher values (see the right 
panel of Fig. [4|). Since we observe multiple recon- 
nection events happening at the same time (compare 
the right panel of Fig. [4] to the Sweet-Parker case in 
Fig. E]), the total reconnection rate should be signif- 
icantly enhanced. 

4.2.1. Dependence on the Turbulent Power 

We run several models with varying power of tur- 
bulence. All other parameters were kept the same. 
This allowed us to estimate the dependence of the re- 
connection rate on the power of injected turbulence. 

Figure [5] shows the evolution of reconnection 
speed in models with the turbulent power P^nj vary- 
ing in the range of over one order of magnitude, from 
0.1 to 2.0. The evolution of Vrec reaches stationar- 
ity in a relatively short period of about one Alfven 
time, estimated from the plot. In order to obtain the 
turbulent power dependence, we averaged Vrec over 
a time interval starting from t = 21 and ending at 
time t = 30. In the subplot of Figure [5] we plot the 
dependence of the averaged reconnection speed on 
the strength of turbulence. Diamonds represent the 



Fig. 6. Time changes of the reconnection speed Vrec for 
models with different injection scale Unj. In the subplot 
we show the dependence of Vrec 

on linj . ^ is the current 
sheet thickness, and it is assumed constant for all models. 

averaged reconnection rate in the presence of turbu- 
lence. Pluses represent the reconnection rate during 
the Sweet-Parker process, i.e. without turbulence. 
The error bars correspond to the standard deviation 
of Vrec^ which is a measure of time variation. 

Fitting to the calculated points gives us a depen- 
dency of the reconnection speed Vrec scaling with the 

1 /2 

power of turbulence as ~ P-^j . The power is pro- 
portional to ~ Vi^/t and t ^ L/Vi - Va/Vi, where Vi 
is the amplitude of turbulence at the injection scale 
and Va is the Alfven speed. This gives the relation 
Pinj ^ thus the dependency of the reconnection 
speed Vrec on the amplitude of fluctuation at the in- 
jection scale Vi is Vrec ^ •, which corresponds to 
the LV99 prediction. 

4.2.2. Dependence on the Injection Scale 

Similar studies have been done in order to derive 
the dependence of the reconnection speed Vrec on the 
scale at which we inject turbulence, Unj . Keeping the 
same power of turbulence for all models we inject 
turbulence at several scales, from kf = 3 to kf = 
12. This limited range allows us to obtain desired 
dependence. 

In Figure [6] we present the results obtained in 
this series of models. From the plot we clearly see 
a strong dependence of the reconnection rate on the 
injection scale. The model with the injection at a 
smaller scale kf = 12 reaches smaller values of the 
reconnection rate. In the model with the injection at 
larger scale /c/ = 3, the reconnection is faster than 
the Sweet-Parker reconnection by a factor of almost 
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10 at some moments. After averaging the rates over 
time, we plot its dependence on the injection scale 
in the subplot of Figure [H 

The fitting to the relation of reconnection rate 
on the injection scale gives the dependency of Vrec ^ 

2/3 1/2 

l-^j , which is stronger then predicted ~ l-^^j by LV99. 
In their paper, the authors considered Goldreich- 



Sridhar model of turbulence (jGoldreich fc Sridhar 
[l995) starting at knj. The existence of the inverse 
cascade can modify the effective knj. In addition, 
reconnection can also modify the characteristics of 
turbulence. This aspect requires more study. 



4.2.3. Dependence on Resistivity 

In the global constraint of reconnection rate de- 
rived by LV99 there is no explicit dependency on the 
resistivity. In order to study this, we performed an- 
other set of models in which we changed the uniform 
resistivity rj only. 
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Fig. 8. The evolution of the reconnection speed Vrec for 
models with different anomalous resistivity ria- In the 
subplot we show the dependence of Vrec on ria- Critical 
value of current density jar it is set to 25.0 in all models. 
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Fig. 7. Dependence of the reconnection rate on the uni- 
form resistivity. Black line shows Sweet-Parker reconnec- 
tion rate, and blue line shows the reconnection rate in 
the presence of turbulence. Dashed lines show actual re- 
connection rates obtained from models. Solid lines show 
corrected speeds according to the procedure described in 
the text. 

In Figure [3 we show reconnection rates obtained 
from this set of models. We plot Vrec for five models 
with T] varying from 2 • 10-"^ to 5 • 10"^. The dashed 
lines show relations obtained directly from simula- 
tions. Since our computational box does not change 
from model to model, the magnetic field has non-zero 
gradients at the boundaries for the higher values of 
resistivity. This affects the evolution by reducing the 
reconnection rate. We know from theory that the re- 
connection rate Vsp during the Sweet-Parker stage 



scales as ~ 77^/^. Using this relation, we have cal- 
culated correction coefficients by taking the ratio of 
T]^^'^ /Vsp. Using these coefficients we correct values 
of reconnection rate obtained during the turbulent 
stage. The result is plotted with the solid blue line. 
The relation signifies virtually no dependence of the 
Vrec on T]. The correctness of this approach requires 
more study. 

In addition to the uniform resistivity dependence, 
we have studied the dependence on the anomalous 
effects as well. The results of these studies are pre- 
sented in Figure [8] and show four models with the 
same uniform resistivity ?7n = 5 • 10~^ and the criti- 
cal current density jcrit = 25.0 but with the anoma- 
lous resistivity parameter rja varying between 0.0 and 
2 • 10~^. In the subplot of Figure [8] we plot the de- 
pendence of the reconnection rate on the anomalous 
resistivity parameter r]a. We see that the reconnec- 
tion speed Vrec is insensitive to the value of r]a to 
within the variations of the reconnection rate in each 
model (see the error bars). 

5. DISCUSSION 

In this paper we tested the model of reconnection 
proposed by LV99. Our results show the significant 
influence of the turbulence on the magnetic recon- 
nection rate. The reconnection speed Vrec (Eq. H]) 
shows dependence on the characteristics of the tur- 
bulence only, the rate of energy injection and the 
scale of injection. In particular, there is no explicit 
dependency on the resistivity. 

Numerical testing of the reconnection model in 
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LV99 is far from trivial. The model is intrinsically 
three dimensional. In order to develop a sufficient 
range of turbulent cascade before it reaches the scale 
of current sheet we have to use high resolution sim- 
ulations. Higher resolution, in addition, minimizes 
the role of numerical diffusion. Another problem is 
the adaptation of proper boundary conditions. Our 
model requires open boundaries in order to allow the 
ejection of the reconnected flux. This type of bound- 
ary conditions should allow for free outflow of the 
matter and magnetic field. This property is crucial 
for the global reconnection constraint, since the re- 
connection stops when the outflow of reconnected 
flux is blocked. 

Even though these numerical simulations allow us 
to study reconnection in the presence of turbulence 
for a limited range of magnetic Reynolds numbers (in 
this paper Rm < 10^), the results provide good test- 
ing of the relations derived by LV99. The strong de- 
pendence of Vrec on the injection scale, when scaled 
to the real conditions of the interstellar medium, 
shows dramatic enhancement of the reconnection 
speed, which even in the presence of a magnetic field 
almost perfectly frozen in the medium, allows for 
fast reconnection with the characteristic time com- 
parable to the Alfven time. The LV99 model pre- 

1 /2 

diets this dependence to be as Vrec ^ ^inj • However, 
our numerical testings show a stronger dependence, 

2 /3 

K-ec ^inj • This difference can be explained by the 
presence of the inverse cascade of turbulence in our 
numerical models. 

The advantage of LV99 model is that it is ro- 
bust and fast in any type of fluid, under the assump- 
tion that the fluid is turbulent. Consequences of this 
are dramatic. The reconnection process is not deter- 
mined by Ohmic diffusion, but it is controlled by the 
characteristics of turbulence, like its strength and 
the energy injection scale. Turbulence efficiently re- 
moves the reconnected flux fulfilling the global con- 
straint (Eq. [Tj). Our numerical results conffim the 
LV99 prediction of a strong dependence of the re- 
connection rate on the amplitude of turbulence, i.e. 
K-ec ^ 1 but most importantly, they show that 
the reconnection in the presence of turbulence is not 
sensitive to the magnetic diffusivity of the medium. 

6. SUMMARY 

In this article we investigated the influence of 
turbulence on the reconnection process using nu- 
merical experiments. This work is the first at- 
tempt to test numer i cally the model presented by 
iLazarian fc Vishniac (|l999l ). We analyzed the de- 
pendence of the reconnection process on the two 



main properties of turbulence, namely, its power and 
the injection scale. We also analyzed the role of 
Ohmic resistivity in the weakly turbulent reconnec- 
tion. We found that: 

• Numerical studies of stochastic reconnection are 
finally possible, even though reconnection in nu- 
merical simulations is always fast. 

• Turbulence drastically changes the topology of 
magnetic field near the interface of oppositely 
directed magnetic field lines. These changes 
include the fragmentation of the current sheet 
which favors multiple simultaneous reconnec- 
tion events. 

• The reconnection rate is determined by the 
thickness of the outflow region. For large scale 
turbulence, the reconnection rate depends on 
the amplitude of fluctuations and injection scale 



as Vr( 



and Vrt 



iflj, respectively. 



• Reconnection in the presence of turbulence is 
not sensitive to Ohmic resistivity. The in- 
troduction of the anomalous resistivity does 
not change the rate of reconnection of weakly 
stochastic field either. 
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